Interacting molecular motors: Efficiency and work fluctuations 
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We investigate the model of "reversible ratchet" with interacting particles, introduced by us earlier 
[Europhys. Lett. 84, 50009 (2008)]. We further clarify the effect of efficiency enhancement due to 
interaction and show that it is of energetic origin, rather than a consequence of reduced fluctuations. 
We also show complicated structures emerging in the interaction and density dependence of the 
current and response function. The fluctuation properties of the work and input energy indicate in 
detail the far-from-equilibrium nature of the dynamics. 
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I. INTRODUCTION 



Molecular motors [lHlOlj are subject to intense study 
both from biological and technological point of view. 
They are paradigmatic examples of machines operating 
at nanometer scale. In a cell, motor proteins powered 
by ATP hydrolysis [HhII help move molecules to places 
where they are needed. Motors assist separation of chro- 
mosomes, copying DN A into RNA and perform many 
more functions |15rll7| . Technological applications of 
the underlying mechanisms flourish [Tel ], including e. g. 
Brownian pumps fl9l . |2~0| and quantum tunneling ratch- 
ets [2l|. They provide also an invaluable testing ground 
for fundamental questions of transport phenomena far 
from equilibrium [221 ] . 

Many models of molecular motors appeared in the lit- 
erature since the pioneering work by Ajdari and Prost 
[23} . The basic mechanism is best elucidated in the mod- 
els which rely on the ratchet mechanism [~], [H, l2~~] - [2~6j 
and also bear the name Brownian motors. The basic 
idea can be viewed either as diffusive motion of a parti- 
cle in spatially asymmetric time-dependent potential or 
as chemically driven transitions between a finite number 
of mechano-chemical states. The former view is more in- 
tuitive, but the latter is closer to reality and opens the 
perspective of fitting the underlying transition probabil- 
ities to experimental data. 

More realistic models are rather built on Markov chains 
in the configuration space constructed as product of spa- 
tial and internal (chemical) coordinates [8|, l27H3"j . This 
approach resides perhaps on more solid experimental evi- 
dence, but the absence of explicit potential makes it very 
difficult to assess the energetic efficiency, the question of 
principal importance in this paper. 

Indeed, one of the points of special interest here will be 
the question of the efficiency of molecular motors. Several 
measures of efficiency can be found in literature. We shall 
use the classical thermodynamic definition r\ — W/E- m , 
where W is work performed and E in energy supplied to 
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the system from external source. Alternative measure 
takes into account viscous resistance (36[, thus reflecting 
better the reality, at the cost that the inequality rj < 1 is 
not guaranteed automatically. Yet other methods of mea- 
suring the efficiency involve explicitly the consumption of 
chemical energy [37|, or the magnitude of the stopping 
force [9]. Note, however, that the former work ( 37]) ex- 
plores the interacting motors and the mechanism if gener- 
ating the non-zero current is related to spontaneous sym- 
metry breaking and this it is principally different from the 
non-interacting case studied in Q. Therefore, the direct 
comparison of the efficiency in these two cases is hardly 
possible. We are not aware of any work in which several 
measures of efficiency would be systematically compared 
on the same model. 

The efficiency of canonical Brownian motors realized as 
either flas hing or rocking ratchets was intensely studied 
[27M29I |38M43| . It turns out that the energetic efficiency 
is rather low[5|, [4l[ , while the experimental data on mo- 
tor proteins, e. g. the kinesin [t| [Hj], report high effi- 
ciency, sometimes even estimated to be close to 100 per- 
cent. We are not in a position to judge the quantitative 
precision of these empirical estimates, although it can be 
suspected that the error margin is rather high. However, 
one is lead to a natural conclusion that the usual ratchet 
mechanism with diffusion as principal driving force is not 
an appropriate model for biological motors. 

In idealized case we can distinguish between ratchet 
and power-stroke mechanisms for molecular motors pij , 
the latter relying rather on quasi-deterministic down- 
hill motion in a free-energy landscape which evolves in 
time. Thus, the particles move as if trapped in a trav- 
eling potential wave. This idea was elaborated in a toy 
model of "reversible ratchet" [HI, H^, |46| , showing much 
higher efficiency, close to the biologically relevant figures. 
Of course, arbitrary combinations and mixtures of the 
ratchet and power-stroke mechanisms can be invented 
and indeed, they are believed to be found in reality, e. g. 
in the myosin V motor (see the review |35| and references 
therein). Nevertheless, it is useful to compare these two 
extremes. We should also note that high efficiency was 
characteristic of the models of either interacting [37| or 
non-interacting (471 |48| motors, which do combine the 
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FIG. 1: (Color online) Graph of the time dependence of the 
potential in which the particles move. Full line: V\ (r), dashed 
line: V2(r). The potential at the third site, Vb, is time inde- 
pendent. 



able on-site repulsion between particles is introduced. We 
found in [72| that not too strong interaction leads to in- 
crease of efficiency. This effect was reproduced qualita- 
tively in analytical calculations based on mean-held (MF) 
approximation. Quantitative agreement was reached in 
an improved MF treatment, developed in (73j . Here 
we investigate the model in depth by further numeri- 
cal simulations. Especially, we elucidate the origins of 
the interaction-enhanced efficiency. We shall show that 
it stems from the energy balance rather than from sup- 
pression of fluctuations. At stronger interaction and/or 
higher density, current reversals and oscillations in re- 
sponse function are found. We also calculate the distri- 
bution of input energy and performed work, which is far 
from being Gaussian. 



ratchet and power-stroke mechanisms. 

The second point we shall concentrate on in this work 
will be the mutual repulsive interaction of molecular mo- 
tors. In the cell, the steric (hard-core) repulsion of motor 
proteins influences significantly their behavior. For ex- 
ample, in gene transcription and translation lar ge n um- 
ber of motor proteins move along the same track ]l5l.[l6|. 
forming the so-called "Christmas tree" structures. Thus, 
they show themselves as a physical realization of the well- 
studied asymmetric exclusion process, introduced first in 
the context of biopolymerization [49|, [5(| and later solved 
exactly in great detail, using sophisticated techniques 
[22I [TO "hi . The model was adapted for molecular mo- 
tors, which may attach and detach with defined rates 
[13, EE]. Later, this situation was studied theoretically 
for the cases of kinesin [56|, ribosomes [13], and RNA 
polymerase 15811 . using the procedures developed in traf- 
fic models [59| . The influence of the geometry of the 
compartment in which the interacting motors diffuse af- 
ter detachment from the track was studied e. g. in [60| . 

Interaction of motors brings about even more compli- 
cated collective effects. In the cell, kinesin and dynein 
molecules typically carry the cargo in groups [8j, |9|, [61| , 
resulting in current reversals. Including explicitly the 
ratchet mechanism of driven diffusion of hard-rod par- 
ticles leads to very intricate effects |62-64|, if the par- 
ticle size and the ratchet periodicity are incommensu- 
rate. The collective movement of coupled Brownian mo- 
tors was studied [65-68] and in some cases the coupling 
was found to induce non-zero current and spontaneous 
oscillations even in mirror symmetric potential due to 
dynamical symmetry breaking [37l . |69| . In analogy with 
these works, the motion of a few rigidly bound motors 
was studied [7(| ■ A special case of such interaction is the 
coordination of the two motor heads within single kinesin 
molecule, which leads also to non-trivial effects 34j. Fi- 
nally, let us mention the interaction of the motors with 
the track, studied in "burnt-bridge" models, e. g. in Ref. 
& 

In our previous paper [72| we introduced a model, 
which is a modified version of the "reversible ratchet". 
Spatial coordinate is discretized, as e. g. in [26]. Tun- 



II. REVERSIBLE RATCHET WITH 
INTERACTING PARTICLES 

Our model contains iV particles occupying integer po- 
sitions on a segment of length L, with periodic boundary 
conditions, and evolves in discrete time. The position of 
i-th particle at the instant t is denoted Xi (t) . The parti- 
cles move under the influence of a variable driving force 
with spatial period 3 and temporal period At. The poten- 
tial of this force is V(x, r) = V(x, r + At) = V x mo d 3(1"), 
at site x and time r. The three independent values of 
the potential V a (T), a = 0,f,2 evolve in a four-stroke 
pattern, with Vq(t) = and the remaining two being de- 
layed one with respect to the other by one quarter-period 
t. Thus, we prescribe 



V 1 (T)=V 2 {r + t) = 
V 

V + 2V {1-r/t) 



- V 

- V 



2V(1 - r/t) 



for 
for 
for 
for 




t 
2t 
3t 



< T < 

< T < 

< T < 

< T < 



t 

2t 
M 
At . 



(1) 



We easily recognize the traveling- wave character of this 
potential, corresponding to the power-stroke mechanism 
of the molecular-motor movement. In all the rest of this 
paper, we fix the amplitude of the potential V = 1. The 
time dependence of the potential is illustrated in Fig. [T] 
Besides the driving potential, there is also a uniform 
external force from the useful load F and, most im- 
portantly, the repulsive interaction from other particles. 
We suppose the interaction is on-site only and we tune 
its strength, in order to interpolate between the non- 
interacting and hard-core cases. The j-th particle feels 
the potential from all remaining ones. To formalize it, 

we denote nj(x, r) = YliLi — j)$( x ~ x i( T )) the num- 
ber of particles, except j-th particle, at site x. (We use 
8 (a — b) for Kronecker delta and 8 (a — b) = 1 — 8(a — b).) 
Thus, the j-th particle moves in the potential 



Uj (x, t) = V(x, t) + x F + 
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n,j(x, 
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FIG. 2: (Color online) Spatio-temporal diagrams of the configurations of the motor particles. The width of the sample is 
L = 18, temperature T = 10, quarter-period t = 4, external load F = 0. Each panel corresponds to different combination of 
two parameters, the number of particles TV and the interaction strength g, whose values are indicated at corresponding columns 
and rows. Dots denote presence of exactly one particle at given space and time, the other symbols presence of more particles, 
namely two (O), three (A), four (o), five (□), and six (v)- The diagonal straight lines are guides for the eye, indicating the 
movement of the minima of the potential V(x,t). 



For g — we recover the non-interacting case, while when 
g — > 1 we approach the hard-core interaction of the ex- 
clusion process [22{ . In contrast with the previous work 
[zH, we use here different form of the interaction in or- 
der to see the limit of hard-core repulsion when g — »• 1. 
Although it may cause some small difficulties when com- 
paring the results of [72| with the present ones, the ad- 
vantage lies in the possibility to see the transition from 
non-interacting case to hard-core repulsion on a finite in- 
terval g G [0, 1]. 

The simulation algorithm proceeds as follows. At each 
integer time r we instantly shift the potential according 
to ([1]). Then, we choose N times a particle randomly and 
let it try to make a jump. Therefore, on average every 
particle makes one attempt per one time unit, but the 
probability that a given particle performs actually k at- 
tempts approaches Poisson distribution with unit mean, 
P(k) — l/(efc!), when N is large. For small N there 



is a deviation from the Poisson distribution, which in- 
duces slight finite-size effects, but in [72| we showed that 
they can be neglected for N larger than about 100. Note 
that in each of these N attempts the external potential 
V(x, t) is the same, but the potential Uj(x,r) felt by 
the particle j may change, because the configuration of 
particles rij(x,T) changes after each particle jump. 

In one attempt, the particle is allowed to jump one site 
to the right or left. The probability of the jump x — > y 
of the j-th particle is 

W^ y = ~ (l + eHudv^-u^r^y 1 (3) 

if \x — y\ = 1 and zero if \x — y\ > 1. For convenience, we 
define the temperature T so that (3 = 270/T. 

Let us now specify the main measurables. The simplest 
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FIG. 3: (Color online) Current (panel (a)) and efficiency 
(panel (b)) as a function of interaction strength, for N = 
L = 1200, t = 16, and F = 0.1. The temperature is T = 
(solid line), T = 10 (dashed line), and T = 30 (dotted line). 



quantity of interest is the current 
J(r)=y( Xi (r+l) 



i(r)) 



(4) 



or rather its time average per particle J = 
lim T _ >00 (T TV) -1 J2t'=iJ( t ')- As we are interested in 
the energetics of the motor, we must define the energy 
input and the useful work done by the particle. The lat- 
ter quantity, at time r, is simply w(r) — FJ(t). The 
energy pumped into the motor from outside relates to 
the change of the potential V a (r) while the particles stay 
immobile. Thus, the energy absorbed by the particle i 
between steps r — 1 and r is 



a i (r) = V(x i {T),T)-V{x i (r),T 
and the efficiency, accordingly, 
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Er'=lEi a i(0 



(5) 



(6) 



Later in this paper we shall investigate the distribution 
of the particle shift 

p ( A - T ) = ji E E 5 { x * ( r + At ) - x * ( r ) - Ax ) ( ? ) 
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FIG. 4: (Color online) Current as a function of interaction 
strength, for L = 1200, t = 16, and T = 10. The number 
of particles is N = 360 (density p = 0.3) in the panel (a), 
N = 840 (density p = 0.7) in the panel (b), and N = 1200 
(density p = 1) in the panel (c). Different curves correspond 
to load , from top to bottom, F — (solid), 0.1 (dashed), 0.2 
(dotted), 0.3 (solid), 0.4 (dashed), 0.5 (dotted). 



and also the joint distribution with the input energy 
P(E iD ,Ax) = 

t+At 



^ -\- ' (8) 

r i r' — r + l 

x S (xi(r + At) — Xi(r) — Ax^j 
where TV is the appropriate normalization. Note that 
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FIG. 5: (Color online) Dependence of the current on the load, 
for L = 1200, N = 600, T = 10, t = 16 and interaction g = 
(•), 0.2 (A), 0.3 (□), 0.4 (O), 0.45 (■), 0.5 (O). 



in both distributions there is implicit dependence on the 
time lag At. 



III. ENHANCED EFFICIENCY 

We show in Fig. [5] examples of typical evolutions of 
the particle configurations, for three densities (p = 1/2, 
p = 2/3, and p = 1) and three interaction strengths 
(g = 0, g = 0.5, and g = 0.9). We can see that without 
interaction, particles are bunched together and dragged 
by the traveling wave of the periodic potential. Interac- 
tion smears out this picture, suppresses the current and 
makes at the same time the local particle density more 
uniform. 

The typical dependence of the current and efficiency 
on the interaction strength is shown in Fig. [3J At zero 
temperature, the dependence contains many steep steps 
with multiple maxima and minima. Therefore, for some 
values of the external load F the current changes sign 
several times when the interaction increases. For larger 
temperatures there are still visible traces of this complex 
dependence, although the singularities (sharp steps) are 
smeared out. We also observe that both the current and 
efficiency approaches zero for very strong repulsion (g — > 
1). We shall see later that this feature is special to some 
values of the particle density p = N/L, for example to 
p = 1, which applies to Fig. [31 The generic feature is 
that for interaction above about g ~ 0.6 the current and 
efficiency approach a constant value. 

The most important finding, from the point of view 
of practical use of the motors, is the increase of the ef- 
ficiency when the interaction is switched on but is not 
too strong. For zero temperature we observe multiple 
maxima of the efficiency, which transform into a unique 
maximum at higher temperatures. The effect of efficiency 
enhancement was investigated in detail in our previous 
work [72j . In this paper we return to the origin of this 
effect later, when we shall discuss the energy balance and 



FIG. 6: (Color online) Dependence of the response at zero 
load on particle density, for N — 1200, t = 16, g = 0.9 and 
the temperature T = 30 (solid line), T = 10 (dotted line), 
and T — (dashed line). 

work fluctuations. 

In Fig. 2] we can see three sets of results for the cur- 
rent, differing in the density of particles. Different curves 
in one set correspond to different external load F. Al 
three cases (and also the data shown in Fig. [3j exhibit 
minimum current, i. e. smallest effective driving, at in- 
teractions somewhere around g ~ 0.4 to g ~ 0.5. In 
order to see what is special in this value of the interac- 
tion, we should note that the change of the potential due 
to presence of a single particle g/(l — g) is equal to the 
amplitude of the traveling-wave potential V = 1 just for 
g = 0.5. At this value of the interaction, one particle 
may block, or at least significantly hinder, the movement 
of the remaining particles. 

We can see that for low density, p < 0.5, the asymp- 
totic current for strongly interacting particles, g — > 1, is 
positive at low load and at the same time is sensitive to 
the value of the load. On the other hand, for 0.5 < p < 1 
the asymptotic current at zero load is negative, i. e, 
the interaction induces current reversal. Contrary to the 
previous case, the asymptotic current seems to be ex- 
tremely weakly dependent on the load. The third panel 
shows again that the asymptotic current is zero for unit 
density, independently of the load. 

Complementary information can be read from Fig. [5J 
showing the dependence of the current on the load. We 
can observe, how the current decreases with the interac- 
tion in the full range of observed F. As a consequence, 
also the stopping force, i. e. the value of F for which 
J = 0, decreases with increasing interaction. It is also 
interesting to note the non-linear decrease of the current 
with the load. So, the response function, defined as the 
derivative dJ/dF, depends on F. 

In Fig. [6] we can see how the zero- load response func- 
tion 



depends on the density, in the regime of very strong 
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FIG. 7: (Color online) Dependence of the response on in- 
teraction strength, for L = N — 1200, and t = 16. Lines: 
response at zero load; symbols: response at zero current. The 
temperature is T — 10 (panel (a)), T = 30 (panel (b)). 



but finite repulsion (g — 0.9). Globally, the response is 
stronger at higher temperature, which is due to the fact 
that at low temperature the movement of the particles is 
determined to larger extent by the traveling wave, with 
lesser influence of the external load, provided the load 
is small. An interesting feature is the structure of the 
peaks and the minima seen in Fig. [6] at all temperatures. 
At integer values of the density the response approaches 
zero. The other minima are not so deep and are located 
at densities slightly above the values p = 1/2, p = 3/2 
etc. Interestingly, the maxima are found at densities very 
close to the fractions p = 1/3, p — 2/3, p — 4/3, and 
p = 5/3. 

As we already said, the response depends on the load, 
so we must distinguish from x F =o a * l eas t one more re- 
sponse function, defined at zero current 



lim ^ 

J^o d J / 



(10) 



We can compare these two quantities in Fig. [7] The 
difference between Xf=o an d X J=0 is especially marked 
for low interaction strength, while at about g ~ 0.3 they 
come close to each other and at g ~ 0.5 the two become 
nearly indistinguishable. The source of this behavior can 
be understood looking at Fig. [SJ Without interaction, 
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FIG. 8: (Color online) Dependence of the current on particle 
density, for N = 1200, T = 10, t = 16, F = and interaction 
g = 0.9 (solid line), g = 0.52 (dashed line), g = 0.49 (dotted 
line), and g — 0.3 (dash-dotted line). 



the dependence of the current on the external load is 
markedly non-linear, so that the derivative at F = and 
J = differ. Increasing the interaction, the non-linearity 
weakens and at g ~ 0.5 we observe nearly linear depen- 
dence, resulting in nearly equal values of the derivative 
at F = and J = 0. Note that the density is p = 1 in 
Fig. [7] and both response functions approach zero when 
g — > 1, in accordance with the results shown in Fig. [6l 

It is interesting that the dependence on g exhibits sev- 
eral peaks. The last (and highest) one is located close to 
g = 1/2 and has nearly the same shape both in x F=0 an d 
Xj=o- However, at lower g the peaks in the two response 
functions are much different. We already mentioned that 
the interaction g = 1/2 is special, as the change in po- 
tential due to presence of a single particle just equals 
the amplitude of the periodic potential V(x). Also the 
second highest peak in x F=0 seems to be located at a 
special value of the interaction strength, namely close to 
g = 1/3. We can also see a small peak close to g = 1/4. 
We believe these special values are due to special blocking 
configurations of particles, which enhance the sensitivity 
of the system to the presence of the external load. In- 
deed, g=l/3 means that two particles on the same site 
contribute as much as the amplitude of V(x), at g = 1/4 
the same holds for three particles at a site. 

To avoid confusion, we do not claim that the config- 
urations of one, two, three, etc. particles are more (or 
less) frequent at certain values of g. What we claim is 
the following. These configurations happen time to time. 
When they do happen, and if g has special values, they 
cause large sensitivity to the value of the load. For other 
values of g, the sensitivity to the load is weaker, whatever 
configuration of particles occurs. 

We also looked at the density dependence of the cur- 
rent at high density and strong interaction. The results 
are summarized in Fig. [5] For the strongest interaction 
investigated, g = 0.9, the curve J(p) has a very pecu- 
liar zig-zag shape, with zeros at p = m/2, maxima at 
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FIG. 9: (Color online) The weight of forward (positive quanti- 
ties) and minus weight of backward (negative quantities) steps 
performed at instants r' = r mod (At) within the period. The 
steps originate at points x' = x mod 3, where x' = corre- 
sponds to dotted line, x' = 1 to solid line, and x' = 2 to 
dashed line. The five panels have the following parameters, 
starting from the top. (a): TV = 400, g = 0, F = 0. (b): 
N = 400, g = 0.9, F = 0. (c): TV = 600, g = 0.9, F = 0. (d): 
TV = 800, 3 = 0.9, F = 0. (e): TV = 400, g = 0.9, F = 0.3. 
All five: L = 1200, T = 10, t = 16. 



p = to + 1/3, and minima at p = m — 1/3, for positive 
integer to. When the interaction is weakened, the sharp 
cusps become mild waves, until the structure of maxima 
and minima vanishes at about g = 0.5. For smaller g, the 
current is a monotonously decreasing function of density. 

Note that the motor with hard-core repulsion under- 
goes a current reversal at a density within the interval 
p G (0, 1). This is in sharp contrast with the asymmet- 
ric exclusion process, where the current is proportional 
to p(l — p). The reason for this difference lies in rather 
different way the particles are driven. In ASEP, there 
is constant and homogeneous drift, only hindered by the 
exclusion principle. In our model, the driving originates 
from the time dependence of the potential, therefore, it 
is also space- and time-dependent. The orientation of 
the current depends on precise timing of the potential 
changes at different places. The interaction changes the 
potential a particle feels and the current is susceptible to 
the details of the potential, so there is no guarantee that 
the orientation of the current will be the same with inter- 



action as it was without interaction. Indeed, the current 
reversal phenomenon demonstrates that the orientation 
of the current does change due to the interaction. Note 
also that the current reversal was observed (qualitatively 
correctly) also in the approximate mean-field calculation 
& 

Some insight into the current reversal phenomenon can 
be gained from the statistics of forward and backward 
steps at different places and different times within the pe- 
riod. We define the measured weight P(Ax = ±1; x', r') 
as the average number of particles which jump forward 
("+" sign) and backward ("— " sign) from site x at time 
r, where x' = x mod 3 and t' = r mod (4t). Note 
that it is not a probability, because it is not normalized 
to unity. We can see a typical example in Fig. |H1 We 
can see that without interaction the particles alternately 
prefer to jump forward from sites x' = 0, 1, and 2. The 
backward jumps are rare. This behavior is independent 
of the particle density by definition. If we add strong 
repulsion, g = 0.9, the picture differs substantially in the 
low and high density regime. For p = 1/3 the statistics 
of forward jumps differs little from the non-interacting 
case, and the frequency of backward jumps is increased, 
but remains low. At half filling, p = 1/2, the parti- 
cles jump alternately forward and backward, at different 
times, so that the total effect is zero current, as seen al- 
ready in Fig. [8l When the density is further increased 
to p = 2/3, the statistics is nearly a mirror image of the 
case p — 1/3. The particles preferably jump backward 
at specific places and times, and the forward jumps are 
rare. For comparison, we show in the last panel of Fig. 
IHl how the statistics is influenced by non-zero external 
load. The time dependence looks similar, but weight of 
forward jumps is suppressed and the weight of backward 
ones is enhanced. 

As the probability of the jumps reflects the local poten- 
tial, and therefore the local instantaneous configuration 
of particles, through the formula ([3]), the statistics of the 
jumps shown in Fig. [5] tells us, what is, on average, the 
local neighborhood of a particle at positions x 1 and times 
t' . Change in the shape of the jump statistics reflects the 
reorganization of the local particle configurations due to 
repulsive interaction. We can clearly see that the reor- 
ganization of the particles can be so dramatic that the 
current changes sign. 

For comparison, we show also the statistics of jumps in 
the presence of non-zero external load. The suppression 
of positive and enhancement of negative peaks is visible, 
as expected. 

Similar analysis can also make more clear the argument 
stated before, that the peaks in the response function at 
special values of g are related to the enhanced sensitivity 
of certain configurations of particles to external load. For 
example, for g — 0.5 such sensitive situation occurs when 
a particle tries to hop to a site where there is already a 
single particle. To support this view we plot a similar 
statistics as in Fig. |H1 but for the difference in the count 
for force F = 0.01 and opposite F = —0.01, on condition 
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FIG. 10: (Color online) Sensitivity of the weight of forward 
(a) and backward (b) steps to the change in external load. In 
this statistics, steps originate at site x, where = x mod 3 
and particles go to the site already occupied by exactly one 
particle. Solid line corresponds to interaction g — 0.5, dashed 
line g — 0.4, dotted line g = 0. Other parameters are N = 
L = 1200, T = 10, t = 16. As for the external load F, see the 
definition of the plotted quantities in the main text. 



that the site to which the particle is moving, already 
contains exactly one other particle. We can write that 
quantity as 

AP(Ax = ±l\n x±1 = 1) = 

= P(Ax = ±lK ±1 = l)| F= _ 001 - (11) 
-P(Ax = ±l\n x±1 = l)\ F=Q01 

where x is the original position of the particle, x ± 1 the 
position after the move, n x ±i number of other particles 
at the site where the particle is about to move. We plot 
an example of this statistics in Fig. [lOj We compare 
the situation at interaction g = 0, g = 0.4, and g = 0.5. 
We can see that the case g — 0.5 is indeed special. The 
sensitivity to the external load is larger. Moreover, the 
difference in statistics has the same sign for almost all in- 
stants within the time period (positive for forward moves, 
negative for backward ones), while both for g = and 
g = 0.4 there are positive as well as negative differences. 
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FIG. 11: (Color online) Current (O) , efficiency (A), average 
energy input (O) and its positive (+) and negative (x) parts, 
all relative to the value at g = 0. The other parameters are 
N = L= 1200, T = 10, t = 16, and F = 0.1. 



IV. WORK FLUCTUATIONS 

To understand better the effect of enhanced efficiency 
due to interaction, we shall look at the energy balance 
and fluctuations. First, we compare the values of cur- 
rent, efficiency and average input energy (-Em) relative 
to their values at zero interaction, denoted Jo, 770, and 
(£'in)o: respectively. More precisely, we plot the typical 
interaction dependence of the quantities J = J /Jo — 1, 

rj = V / m - 1 and = (E in )/(E in ) - 1 in Fig. [JU 

We can clearly see that the initial increase of efficiency 
for small g is accompanied by nearly no change in the 
current, while the input energy decreases. Therefore, the 
enhanced efficiency is due to lower energy input, while the 
output (the work) remains nearly unchanged. When the 
interaction strength increases further, the current starts 
decreasing as well and so does the work, which is pro- 
portional to J. This effect finally outweighs the lower 
energy input and the efficiency decreases again. This is 
the source of the maximum in the efficiency at moderate 
values of the interaction. 

We can get a bit more detailed information if we split 
the input energy into its positive and negative parts. Re- 
call, that according to the definition ((5J) the input energy 
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FIG. 12: (Color online) Dependence of the average energy 
input gained by one particle on the shift of the particle. The 
quarter-period is t = 16 and the time lag At = lOOi (i. e. 
25 periods). The temperature is T = 10 (O and A), T = 30 
(□ and O). The interaction is g = (O and O), g = 0.1 
(□ and A). The remaining parameters are N = L = 1200, 
F = 0.1. The horizontal axis is normalized by the number of 
time periods, which is At /it, therefore it expresses the shift 
per period. 
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FIG. 13: (Color online) Probability distribution of the input 
energy gained by a single particle at fixed value of the shift 
of this particle, within the time lag At — 400 1 (i. e. 100 
periods). The fixed shift is Ax = 282 (O), Ax = 283 (□), 
and Ax = 284 (A). The other parameters are N = L = 1200, 
T = 10, t = 16, g = 0.1, and F = 0.1. The horizontal axis is 
normalized by the number of time periods, which is At /At, 
therefore it expresses the energy input per period. 



is E- m = X)t'=t+i a i( T ')- We separate the contributions 
from times when aj(r') is positive from those when it 
is negative. The former correspond to the shift of the 
potential V(x,t) upward, that latter to its downward 

move. So, E in± = J2l t=J +l ±a i (r / )0( ± ai(r')) , where 
6{x) is the Heaviside function. With this definition we 
have E- m — E- ln + — E- lu -. We then define the contribu- 
tions from positive and negative potential moves to the 
quantity (E in ) as (E in ) + = ({E- m+ )-(E in+ } )/(E in ) and 

(E in )_ = ((Ejn-)o - (£ in -))/(£m)o where, as above, the 
subscript denotes the quantities computed at g = 0. 
Hence (E in ) = (E in ) + + {E in }_. We show the depen- 
dence of (Ei n ) ± again in Fig. [TT] We can see that both 
positive and negative parts contribute to the decrease of 
the input energy. The contribution of the positive part is 
larger in the most interesting region of moderate g, where 
the efficiency grows with interaction, while for larger g 
the decrease of the negative part becomes more impor- 
tant. This leads to the following explanation of the effect 
of increased efficiency. 

At not too high temperature, the particles are chiefly 
driven by the traveling wave of the periodic potential. 
This is the power-stroke mechanism of the molecular mo- 
tor. When the interaction is switched on, but remains 
small, the particles move in an effective potential which 
differ little from the original traveling wave. So, the cur- 
rent remains nearly the same, as testified in Fig. 1111 
while the input energy is lowered, as is also seen in Fig. 
ITT1 This lowering could be understood as follows. On 
the other hand, the repulsion affects the distribution of 
particles within the period of the potential V(x). The 
minima become shallower, therefore the particles are less 



concentrated at them. But it is the minimum of the 
potential which is shifted above in the time evolution 
of the potential, so it is the particle at the minimum 
that acquires the energy from the source of the driving. 
Less particles at the minimum equals less input of energy, 
more precisely lowering of the positive part of the input 
energy. Conversely, the particles pushed off the instanta- 
neous minima of the potential are found at the maxima 
of the potential. But these particles suffer lowering of 
the time-dependent potential, i. e. returning the energy 
back to the external source, therefore lowering also the 
negative part of the input energy. These two effects, i. 
e. unchanged current and lowered energy input, are the 
explanation of the increased efficiency. Of course, more 
subtle effects are also at work here. Especially, also the 
negative part of the energy input contributes. More im- 
portantly, if the interaction is strong enough, it changes 
the potential the particles move in to such extent that 
the current diminishes. At very small temperature, the 
current is sensitive to tiny changes in the shape of the 
potential and small changes in the interaction strength 
can cause big jumps in the current. We have seen these 
jumps in Fig. [3] 

In addition to the averages, we measured also the full 
joint distribution function of particle shift and input en- 
ergy ([5]). Because the work performed by one particle is 
proportional to its shift, we have in fact the joint distri- 
bution of performed work and input energy. As a first 
piece of information we plot in Fig. the average en- 
ergy input at fixed value of the particle shift, during the 
time interval At. We can observe the already discussed 
fact that interaction decreases the energy input. Here we 
can see that it holds also for most values of the shift, i. 
e. work performed by one single particle, separately. 
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FIG. 14: (Color online) Probability distribution of the shift of 
a particle within the time lag At — 1000 1 (i. e. 250 periods), 
relative to the most probable value Axq. The interaction 
strength is g = (O, with Axo = 718) and g = 0.1 (A, with 
Ax = 718). The other parameters are N = L = 1200, T = 
10, t = 16, and F = 0.1. The horizontal axis is normalized 
by the number of time periods, which is At /At, therefore it 
expresses the shift per period. 



The probability distribution of the energy input at 
fixed shift is shown in Fig. [13J We can see that the 
shape is pretty close to a Gaussian. This is far from 
being true for the distribution of the shift, which is pro- 
portional to the work performed by a single particle, as 
shown in Fig. 1141 The distribution is skewed; when 
we compare the shifts shorter and longer than the most 
probable value, we find that the shorter are significantly 
more probable than the longer ones. This is due to the 
far-from-equilibrium character of the process. We can 
also see that the distribution is composed of two sep- 
arate branches. The first one, with higher probability, 
corresponds to shifts which are multiples of 3, the period 
of the potential. The other shifts have significantly lower 
probability. In fact, it comes as no big surprise, that af- 
ter integer number of time periods the particles like to 
be shifted by integer number of spatial periods. 

The most important finding, however, stems from the 
comparison of the distribution in the case with and with- 
out interaction. In Fig. [T3]we make this comparison for 
such set of parameters where we know that the interact- 
ing case exhibits higher efficiency. By analogy with equi- 
librium statistical physics one might be tempted to guess 
that higher efficiency is accompanied, or even caused, by 
milder fluctuations. The opposite holds in our model of 
the molecular motor. The fluctuations of the work per- 
formed by a single particle are higher in the interacting 
case. Therefore, we conclude that the increase of effi- 
ciency is not accompanied by decrease of fluctuations. 
On the contrary, the study of the energy balance dis- 
cussed above together with the fact of increased fluctua- 
tions shows that the enhancement of efficiency is purely 
an energy effect. 



V. CONCLUSIONS 

Interacting molecular motors moving in the power- 
stroke regime were modeled using a "reversible ratchet" 
model. Tunable on-site repulsive interaction leads to a 
host of intricate phenomena. The most important of 
them is the increase of energetic efficiency for small to 
moderate values of the interaction strength. We traced 
the origin of this effect to energy balance of the process. 
When the interaction is increased from zero, the per- 
formed work remains practically unchanged, while the 
input energy decreases. At the same time, the fluctua- 
tions of the performed work increase. This implies that 
the enhanced efficiency does not originate from the sup- 
pression of fluctuations, contrary to the situation in equi- 
librium heat machines. 

Moreover, we observed that at very low temperature 
the dependence of current as well as efficiency on the in- 
teraction strength is rather complex, composed of many 
upward and downward steps. Hence, the efficiency has 
several, rather than single, local maxima as a function of 
interaction. As for the current, for suitable values of the 
parameters we can observe a sequence of current reversals 
when we increase the interaction strength. This compli- 
cated behavior is due to the interplay between size of 
steps in the external periodic potential, in which a par- 
ticle moves, and the size of additional contributions to 
the potential from the interaction with other particles. 
However, this complicated dependence gradually disap- 
pears when the temperature increases. But the effect of 
current reversal due to interaction remains still visible. 

We also investigated the response function of the cur- 
rent with respect to external load, both for zero load 
and for zero current. We showed that these two response 
functions differ substantially at zero or small interaction, 
but become identical when the interaction is large. We 
also revealed the structure with several peaks for both 
density and interaction dependence of the response func- 
tion. Detailed study of the location of these maxima 
and minima showed that they correspond to specific frac- 
tional values of density and interaction. For example, the 
response is zero if the density is integer number and has 
maximum for densities equal to integer number of thirds, 
except the values which are themselves integers. In the 
interaction dependence, the peaks were found close to in- 
teraction strength equal to one half, one third and one 
fourth. We speculate that these special values are due 
to the fact that in those cases just one, two, and three 
particles on the same site, respectively, contribute to the 
potential by the value exactly equal to the amplitude of 
the external periodic potential. Contrary to the compli- 
cated step structure in the current, the peak structure 
in the response function survives also at higher temper- 
atures. 

The probability distribution of performed work and 
input energy reveals that the interaction leads to the in- 
crease of fluctuations, as we already mentioned. But we 
can see more. First, the distribution of work is far from 
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Gaussian. It is skewed so that the lower particle shift (i. 
e. work performed by an individual particle) relative to 
the maximum is more probable. This is the sign of far- 
from-equilibrium regime of the transport in the molecular 
motor. On the other hand, the input energy is Gaussian- 
distributed, when observed at fixed work. 

There is also a very interesting principal question re- 
lated to large-deviation properties of the fluctuation of 
the performed work. We made some simulations in this 
direction, which show that the work distribution, when 
properly rescaled, converges to a large-deviation func- 
tion. In the last decade, there was a great surg e of ac- 
tivity in the field of Fluctuation Theorems [741484 but 
in our case the problem of applying these results lies in 
the choice of the proper quantity which would be both 
physically meaningful (or at least the physical meaning 
must not be enormously intricate) and satisfy the Fluc- 
tuation Theorem in some of the forms known so far. This 
question remains open. 

Finally, we must also admit several drawbacks of our 
model, which can be lifted only by setting up a com- 
pletely different scheme of particle movement. The first 
point is that the potential changes synchronously at all 
sites. This is unrealistic in biological motors, where each 
molecule has its internal "clock" telling in what phase of 
the chemical cycle the motor finds itself. It could be eas- 



ily possible to simulate an asynchronous version of the 
algorithm. On the other hand, in technological appli- 
cations the synchronicity in the potential changes may 
be built in into the system. The second point concerns 
the tunable interaction used in our model. Motor pro- 
teins interact always by steric hard-core repulsion and 
th effective weak repulsion may occur only as projection 
of real three-dimensional situation onto one-dimensional 
effective model [85|,[86j]. However, there is no principal 
obstacle to simulate three-dimensional case directly, if 
only sufficient computer power is available. Another way 
out is to generalize the asymmetric exclusion process in 
such a way that the maximum number of particles one 
site may accommodate is not one, but two, or three, etc. 
Simulations in this direction are under way. 
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